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Abstract 

In view of ongoing experimental activities to determine the pion polarizabilities, 
we have started to recalculate the available two-loop expressions in the framework 
of chiral perturbation theory, because they have never been checked before. We 
make use of the chiral Lagrangian at order now available, and of improved tech- 
niques to evaluate the two-loop diagrams. Here, we present the result for the neutral 
pions. The cross section for the reaction 77 — ► 7r°7r° agrees with the earlier calcu- 
lation within a fraction of a percent. We present analytic results for the dipole and 
quadrupole polarizabilities, and compare the latter with a recent evaluation from 
data on 77 — ► 7r°7r°. 

PACS: 11.30.Rd; 12.38.Aw; 12.39.Fe; 13.60.Fz 

Key words: Chiral perturbation theory; Two-loop diagrams; Pion polarizabilities; 
Compton-scattering 



1 Introduction 



We consider the process 77 — > 7r 7r° in the framework of chiral perturbation 
theory (ChPT ) [1,2]. The one-loop calculation of the scattering amplitude was 
performed in Refs. [3,4], and the two-loop amplitude was worked out in [5]. 
Because the effective Lagrangian at order p 6 was not available at that time, the 
ultraviolet divergences were evaluated in the MS scheme, then dropped and 
replaced with a corresponding polynomial in the external momenta. The three 
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new counterterms which enter at this order in the low-energy expansion were 
estimated with resonance saturation. Whereas such a procedure is legitimate 
from a technical point of view, it does not make use of the full information pro- 
vided by chiral symmetry. The evaluation of the two-loop amplitude involving 
charged pions was performed later by Burgi [6]. 

Over the last ten years, considerable progress has been made in this field, both 
in theory and experiment. As for theory, the Lagrangian at order p 6 has been 
constructed [7,8], and its divergence structure has been determined [9]. This 
provides an important check on the above calculations: adding the counterterm 
contributions from the p 6 Lagrangian to the MS amplitude evaluated in [5] and 
in [6] must provide a scale independent result. Also in the theory, improved 
techniques to evaluate the two-loop diagrams that occur in these amplitudes 
have been developed [10]. The improvement arises mainly in the evaluation of 
diagrams with four external legs, where the techniques of Ref. [10] allow one 
to extract the ultraviolet divergences by use of simple recursion relations. We 
are now able to present the final result for the two-loop amplitudes in a rather 
compact form (in Refs. [5,6], the result was presented partly in numerical form 
only, because the algebraic expressions were too long to be published). 

Concerning experiment, quadrupole polarizabilities [11] for the neutral pions 
have recently been determined from data on 77 — > 7r°7r° [12]. Further, the 
charged pion polarizabilities (a — (3) n + have been determined at the Mainz 
Microtron MAMI [13], with a result that is at variance with the two-loop 
calculation presented in [6]. Last but not least, there is an ongoing experiment 
by the COMPASS collaboration at CERN to measure the charged pion and 
kaon polarizabilities [14,15]. 

In view of these developments, and because the two-loop expressions for the 
polarizabilities had never been checked, we decided to recalculate these ampli- 
tudes, using the improved techniques of Ref. [10] to evaluate the integrals, and 
invoking the chiral Lagrangian at order p 6 [8,9]. As the calculation in the case 
of neutral pions involves considerably less diagrams, and because the Fortran 
code for these amplitudes is still available to us for checks, we have decided to 
start the program with a re-evaluation of these amplitudes. This is the main 
purpose of the present work. The evaluation of the corresponding expressions 
for the charged pions and for the kaons is underway and will be presented 
elsewhere [16]. 

The article is organized as follows. Section 2 contains the necessary kinematics 
of the process 77 — > 7r°7r°. To make the article self contained, we summarize in 
Section 3 the necessary ingredients from the effective Lagrangian framework. 
In Section 4, we display the Feynman diagrams and discuss their evaluation. 
Section 5 contains a concise representation of the two Lorentz invariant am- 
plitudes that describe the scattering matrix element. In Section 6, we compare 
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the present work with the previous calculation [5] , while Section 7 contains ex- 
plicit expressions for the dipole and quadrupole polarizabilities valid at next- 
to-next-to-leading order in the chiral expansion, together with a numerical 
analysis and a comparison with an evaluation from data on 77 — > 7r°7r° [12]. 
The summary and an outlook are given in Section 8. Finally, several technical 
aspects of the calculation are relegated to the Appendices. 



2 Kinematics 

The matrix element for the reaction 

7(gi)7fe)- A Pl )7r> 2 ) (2.1) 

is given by 

(vrV) out I 7 (?i)7(fe) W = < (2tt) V 4 ) (P f - P t ) T N , (2.2) 

with 

1 — e fc l fc 2 v V-v 1 

V^ = ij dxe^ x+ ™\Tt\ Vl )ix\ V2 ) out I Tj,(x) 3l/ (y) | 0). (2.3) 

Here is the electromagnetic current, and a = e 2 /47r ~ 1/137. We consider 
real photons, q 2 = 0, with q • qi = 0. The decomposition of the correlator V^ u 
into Lorentz invariant amplitudes reads 

Vp, = A(s, t, u)T 1>lv + B(s, t, u)T 2ilu + C(s, t, w)T 3 ^ + D(s, t, u)T 4fM/ , 
T± = - s g^ — q\ v q 2lx , 

T 2fMU = 2s A^A U - v 2 g^ ~2v(q lu A^ - q 2fl A u ) , 

Ti^ = s(q lfl A u - q 2u A^) -v{qi^qi u + q 2 ^q 2u ) , 
A M = ( Pl -p 2 ) M , (2.4) 
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where 



s = (<?i + g2) 2 , t=(pi-qi) 2 , u = (p 2 -qi) 2 , v = t-u (2.5) 

are the standard Mandelstam variables. The tensor V^ u satisfies the Ward 
identities 

<fiV lu , = <£V IM , = 0. (2.6) 

The amplitudes A and B are analytic functions of the variables s, t and u, 

symmetric under crossing (t, u) — > (u,t). The amplitudes C and D do not 
contribute to the process considered here, because q • = 0. 

It is useful to introduce in addition the helicity amplitudes 



H ++ = A + 2(4M%-s)B, iJ + _ = 8(M ' tU) B. (2.7) 

s 

The helicity components -£/++ and H + - correspond to photon helicity dif- 
ferences A = 0,2, respectively. With our normalization of states (P1IP2) = 
2 (2%) 3 p® 5^(pi — p 2 ), the differential cross section for unpolarized photons 
in the centre-of-mass system is 



— = — 0(a) lf(M) , #(M) = |# ++ | 2 + I# + -| 2 , (2.8) 



with f3(s) = y 1 — 4M^/s. The relation between the helicity amplitudes M + ± 
in Ref. [12] and the amplitudes used here is 



M++(s, t) = 27raH ++ (s, t) , M+_(s, *) = 16naB(s, t) . (2.9) 

The physical regions for the reactions 77 — > 7r°7r° and 771° — > 771° are displayed 
in Fig. 1, where we also indicate with dashed lines the nearest singularities 
in the amplitudes A and B. These singularities are generated by two-pion 
intermediate states in the s, t and u channel. 



3 The effective Lagrangian and its low energy constants 



The effective Lagrangian consists of a string of terms. Here, we consider QCD 
with two flavours, in the isospin symmetry limit m u = m d = m. At next-to- 
next-to-leading order (NNLO), one has 
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Fig. 1. Mandelstam plane with three related physical regions, s-channel: 77 — > irn, 
t- and u-channel: jtt — ► 77r. We denote the threshold for 77 — > 7T7r (77r — ► 77r) by 
A (P). The dashed lines at s,t,u = 4M 2 indicate the presence of branch-points in 
the amplitude, generated by two-pion intermediate states. 

C cS = C 2 + C 4 + C 6 . (3.1) 
The subscripts refer to the chiral order. The expression for C 2 is 

F 2 

C 2 = —(DJJDorf + M 2 (U + 17*)) , 
D,U = d„U - i{QU - UQ)A, , Q = |diag(l, -1) , (3.2) 

where e is the electric charge, and denotes the electromagnetic field. The 
quantity F denotes the pion decay constant in the chiral limit, and M 2 is 
the leading term in the quark mass expansion of the pion (mass) 2 , M 2 = 
M 2 (l + 0(m)). Further, the brackets (. . .) denote a trace in flavour space. In 
Eq. (3.2), we have retained only the terms relevant for the present application, 
i.e., we have dropped additional external fields. We choose the unitary 2x2 
matrix U in the form 



7T 



U = a + in/F , a + — = l 2x2 , tt 



1 7T° V27V- 



V2 



(3.3) 



7T — 7T 



The Lagrangian at NLO has the structure [1] 
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U = E^ = jiD^U D"rt) 2 + • • • , (3.4) 
i=i 4 

where denote low energy couplings (LECs), not fixed by chiral symmetry. 
At NNLO, one has [8,9] 
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£ 6 = 5>Pi. (3.5) 

i=i 

For the explicit expressions of the polynomials Pi, we refer the reader to 
Refs. [8,9]. The vertices relevant for 77 — > 7r°7r° involve 1±,...,Iq from £ 4 
and c 29 , . . . , C34 from £ 6 . 

The couplings ^ and q absorb the divergences at order p 4 and p 6 , respectively, 



= c)*" 4 d) + 7i A}, 
d = d) - 7? } A 2 - (7f } + 7 f } (ix, d)) A} , 

A = Ij ^^ > lnc=-i{ln4,r + Hl) + l}. (3-6) 

The physical couplings are /[(//, 4) and c[(/i, 4), denoted by /[,c[ in the fol- 
lowing. The coefficients % are given in [1], and t^ 1 ' 2 ' 1 ^ are tabulated in [9]. In 
order to compare the present calculation with the result of [5], we will use the 
scale independent quantities k introduced in [1], 



where the chiral logarithm is / = ln(M^//z 2 ). We will use [17] 



h = -0.4 ± 0.6 , l 2 = 4.3 ± 0.1 , l 3 = 2.9 ± 2.4 , l 4 = AA ± 0.2 , (3.8) 
and [18] 

J A = i 6 -J 5 = 3.0 ± 0.3 . (3.9) 
The constants c\ occur in the combinations 



^ = 40967r 4 (-c r 29 -cS + cy , 
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a r 2 = 256vr 4 (8 c r 29 + 8 c r 30 + c r 31 + c r 32 + 2 c r 33 ) , 



128tt 4 (c r 31 + c r 32 + 2 c£ 




(3.10) 



Their values have been estimated by resonance exchange e.g. in Ref. [5], see 
also [19], where c 34 has been determined from a chiral sum rule. For the present 
application, we simply take the values obtained in [5], 



In the numerical evaluations discussed later on, we use /z = M p . As mentioned 
already in Ref. [5], varying this scale between 500 MeV and 1 GeV leads to 
a negligible change of e.g. the cross section for the reaction 77 — > 7r°7r° below 
400 MeV. Finally, we will use F n = 92.4 MeV [20] (see [21] for a recent update 
of this value), and M n = 135 MeV. 



4 Evaluation of the diagrams 

The lowest-order contributions are the one-loop diagrams displayed in Fig. 2. 
They have been evaluated for the first time in Refs. [3,4], where it was noticed 
that the sum of these two amplitudes is ultraviolet finite, because there are 
no contributions from the effective Lagrangian at order p 4 at this order. The 
two-loop diagrams are displayed in the Figs. 3,5 and 6. The two-loop diagrams 
in Fig. 3 may be generated according to the scheme indicated in Fig. 4, where 
the shadowed blob denotes the ci-dimensional elastic 7T7r-scattering amplitude 
at one-loop accuracy, with two pions off-shell. As is discussed in Appendix B, 
the one-loop integrals in the nn amplitude may be represented in a dispersive 
manner. This allows one to reduce the two-loop integrals in Fig. 3 to the one- 
loop ones, where one has to perform at the end an integration over a dispersive 
parameter. 

Two further diagrams are displayed in Fig. 5. The first one - called "acnode" 
in the literature - may again be evaluated by use of a dispersion relation, see 
Appendix B. The second one is trivial to evaluate, because it is a product of 
one-loop diagrams. The remaining diagrams at order p 6 are shown in Fig. 6. 

The evaluation of the diagrams was done in the following manner. 

(1) We have performed the integration over the loop momenta in the ci- 
dimensional regularization scheme, in particular using the procedure de- 
scribed in Ref. [10] and in Appendix B, and invoking FORM [22]. 



a\{M p ) + 86 r (M p ) = -14±5, 
a r 2 {M p ) -2b r {M p ) = 7±3, 

b r (M p ) = 3 ± 1 ; M p = 770 MeV. 



(3.11) 
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Fig. 2. The one- loop diagrams. 




(1) (2) (3) (4) 




(5) (6) (7) (8) 




(9) (10) (11) (12) 




(13) (14) (15) (16) 




(17) (18) (19) 



Fig. 3. A set of two-loop diagrams generated by £2 and one-loop diagrams generated 
by £4. 

(2) We have then checked numerically that the amplitude satisfies the Ward 
identities (2.6) in d dimensions, in the unphysical region, where the am- 
plitudes are real. 

(3) We have verified that the counterterms from the Lagrangian [9] re- 
move all ultraviolet divergences, which is a very non-trivial check on our 
calculation. 

(4) We have checked that the (ultra-violet finite) amplitude so obtained is 
scale independent. 

(5) Finally, we have numerically verified that the three lowest partial waves of 
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<?1 



+ 



S d d l 




/ 



qi+l 



qi — I 




Fig. 4. Construction scheme for the diagrams in Fig. 3. 





(1) (2) 

Fig. 5. Acnode and butterfly diagrams. 



\ 



, 2 'N/V^/V 1, 2 ^^^\/V/N/\, 




\ 



29 - 34 



(1) 



(2) 



(3) 



(-1) 



Fig. 6. The remaining diagrams at order p e : one-loop graphs generated by £4, and 
counterterm contributions from L%. 



the helicity non-flip amplitude H ++ carry the proper one-loop tttt phase, 
in conformity with unitarity. 

We note that the steps (3) and (4) could not be performed in Refs. [5,6], 
because the counterterms at order p 6 were not yet available [9] . 
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5 The two-loop amplitudes 



We give the expression for the amplitudes A and B by using the same notation 
as in [5] , and refer the reader to Appendix C of this reference for the one- loop 

integrals J{s), J {s),G{s),G (s) and H{s) that occur. We have 



A = ^^-(s-M 2 J + U A + P A + 0(E*), 

S TV 

B = U B + P B + 0(E 2 ). 



(5.1) 



The quantity G n (s) stands for G(s), evaluated with the physical pion mass. 
The unitary parts Ua(b) contain s,t and ^-channel cuts, and Pa(b) are linear 
polynomials in s. We find 



U A 



{k - 5/6) 



144tt 2 sF4 
+A A (s,t,u) , 



G(s) [(s 2 -Mt)J{s) + C{sJ t ) 
s-AMl) {h(s)+ 4 



A 



(s - Ml) J{s) 



24vr 2 F4 
sG(s)+2M 2 (G(s)-3 j(s)) 



d 2 



C(s,h) = {2 («i — |) (*-2M 2 ) 2 + I (/ 2 -£) (4s 2 - 8sM 2 + 16 M*) 



-3M^r 3 + 12M 2 (s-M 2 )r 4 - 12sM 2 + 15M^} , 



4 = ^(3 cos 2 ^-l), 



(5.2) 



(l 2 - 5/6) H(s) . , 

The expressions for A^(b) are displayed in the Appendices C and D. 
The polynomial parts are 



(5.3) 



1 



a-i M 2 + a 2 s 



A (16tt 2 F2) 2 

a x = a\ + ^^Al 2 + 1 (%T 2 + 12T&-- 

1 f / - 4\ 5 - 6971 

° 2 = ° 2 " 18 V + ^ [ 2l2 + 12[a " 3) " 3 ' 2 + 12h + 144/ ' 



20 7 , rtT 1101 
y / 2 + 12/ A + -}, 



(5.4) 
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Pb 



b 



(16 TT^) 



2 ' 



b 




393" 
144. 



(5.5) 



The constants a\,d' 2 and b r are displayed in terms of the LECs at order p 6 
in Eq. (3.10). Using the fact that the bare couplings q displayed in Eq. (3.6) 
are scale independent, one indeed finds that the above expressions for the 
amplitudes A, B are scale independent as well. 



6 Comparison with the previous calculation 

We can now compare the amplitudes A, B with the earlier calculation, pre- 
sented in Section 7 of Ref. [5]. In that reference, the amplitudes were evaluated 
with a different techniques. Furthermore, the Lagrangian was not available 
in those days, and an important ingredient to check the final result was, there- 
fore, missing. We can make the following observations. 

(1) The amplitudes A and B consist of a part with explicit analytic expres- 
sions, and additional terms Aa,b, that are given in the Appendices C and 
D of the present work in the form of integrals over Feynman parameters. 
These latter terms were given only in numerical form in [5]. 

(2) The explicit analytic expressions agree with the previous calculation, ex- 
cept for the coefficient of the single logarithm in 02 in Eq. (5.4). The 
factor 2/3 in [5] is replaced by - 4/3 here. As the present amplitude is 
scale independent, we conclude that it is the result Eq. (5.4) which is 
correct 1 . This mistake does not affect the algebraic expressions for the 
polarizabilities discussed below, for which we fully agree with Ref. [5]. 

(3) We can compare the quantities Aa,b in numerical form only. For this 
purpose, we have made two checks. First, we have evaluated the cross 
section for the reaction 77 — > 7r°7r° below a centre-of-mass energy of 400 
MeV, using the same values for the LECs as in [5]. It agrees with the 
previous one within a fraction of a percent - the difference would not 
be visible in Fig. 5 of Ref. [5], and we do not, therefore, reproduce that 
plot here. Second, we have re-evaluated the two-loop contributions to 
the polarizabilities presented in column 4 of Table 3 in [5] . The numbers 
(0.17, - 0.31) in the old calculation become (0.17, - 0.30) here. 

1 Burgi [23] provides in his thesis work the isospin / = 0, 2 amplitudes, and the 
one for the charged pions. Subtracting the latter from the former reveals that his 
calculation agrees with the statement just made. 
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(4) To summarize, we confirm the previous result up to the coefficient in one 
of the chiral logarithms, and up to minute changes in the numerical values 
of Aa,b- Numerically, the results in [5] are not affected in any significant 
manner by these modifications, whose effect is by far smaller than the 
uncertainties generated by the (not precisely known) values of the low 
energy constants. 



7 Pion polarizabilities: dipole and quadrupole 



The dipole and quadrupole polarizabilities are defined [11,12] through the ex- 
pansion of the helicity amplitudes at fixed t = M 2 , 



— H +T (s, t = Ml) = (a, + ft)*,, + —(a 2 ± 0O.o + 0(s 2 ) . (7.1) 

Because we have at our disposal the helicity amplitudes at two-loop order, we 
can work out the polarizabilities to the same accuracy. It turns out that all 
relevant integrals can be performed in closed form. We discuss the results in 
the remaining part of this Section. 



7.1 Chiral expansion 

Using the same notation as in [5], we find for the dipole polarizabilities 
16tt 2 F2M w \ 16tt 2 F 2 

with 



(«i ± = TT-^TTT U± + + 0{Mt)\ , (7.2) 



Ci + = 0,C!_ = -1/3, 

d 1+ = 8 b r - — (144 1(1 + 2 1 2 ) + 96 / + 288 1 2 + 113 + A+) , 
648 

d x _ = a[ + 8 V + J- (144 / (3 l A - 1) + 36 (8 h - 3 Z 3 - 12 f 4 + 12 l A ) 
648 ^ 

+43 + A_) , 

A+ = 13643 - 1395 tt 2 , A_ = -3559 + 351 tt 2 . (7.3) 

We have split off the numbers 113 and 43, respectively, to illustrate that these 
expressions completely agree with the ones displayed in Eq. (8.14) of Ref. [5], 
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Table 1 

The dipole and quadrupole polarizabilities in units of fm 3 and 10~ 4 fm 5 , re- 
spectively. The symbol [*] refers to the present work. The slight difference in the 
value of («i + Pi) n o with the one reported in [5] is due to the updated values of l 2 
and of the pion decay constant F w used here. 





ChPT 


dispersion relations 


("1 - 


-1.9 ±0.2 [5] 


-1.6 ±2.2 [24] 


("1 + 0l)n° 


1.1 ±0.3 [5] 


0.98 ±0.03 [24] 
1.00 ±0.05 [25] 


(« 2 - 


37.6 ±3.3 [*] 


39.70 ± 0.02 [12] 


(«2 + 02)ir° 


0.04 [*] 


-0.181 ±0.004 [12] 



where, as already mentioned, no explicit expressions for the remainders A + _ 
were worked out. For the quadrupole polarizabilities, we obtain 



(a 2 ± P 2 ) n o = 



a 



16 tt 2 Fl Ml 



c 2 ± + 



Mid 



2± 



16 7T 2 Fl 



+ 0(M 4 ) , 



(7.4) 



with 



c 2+ = 0, c 2 _ = 156/45, 
d 2 + ~- 



5009 13453 tt 2 16Z 2 
+ + 



27 



720 



45 



= 12 a r 2 - 24 W + — (l280 1 (1 - 6 l A ) + 19216 - 181 1 tt 2 ) 



960 

4 (52 h + 5 T 2 + 3 k ~ 78 T 4 + 105 l A ) 
45 



(7.5) 



(7.6) 



7.2 Numerical results 



For numerical evaluations of the polarizabilities we use the values of the LECs 
given in Section 3. The results are displayed in Table 1, where we also quote 
the results from dispersive calculations. The following comments are in order. 

(1) The slight difference in the value of (aq + /? 1 ) 7r o with the one reported in 
[5] is due to the updated values of l 2 and of the pion decay constant F n 
used here. 
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2 
1 


CO 

'q-2 



-3 
-4 



" 5 -2 -10 12 

sM" 2 

Fig. 7. The helicity amplitudes j^H ++ (s,t = Ml) (solid line), j^H + _(s,t = Ml) 
(dashed line), plotted as a function of s, at t = M%. Compare with the definition of 
the polarizabilities in Eq. (7.1) The dotted lines are displayed to guide the eye. 

(2) Our results for the dipole polarizabilities as well as for the quadrupole po- 
larizabilities («2— ^2)77° agree with the results of the recent investigations 
performed in [25,24] and [12] within the uncertainties quoted. 

(3) The prediction for the quadrupole polarizability (02 + ^2)1^ is positive, in 
contrast to the result reported in Ref. [12]. The ChPT expression contains 
as the only LEC I2, known rather accurately from tttc scattering [17]. We 
come back to this point in the following subsection, where we also discuss 
the uncertainties quoted in the Table for the ChPT calculation. 

(4) We plot the helicity amplitudes in Fig. 7. It illustrates the fact that the 
helicity flip amplitude H + _ is quite flat at this order, in contrast to the 
non-flip amplitude H ++ , see the values of the quadrupole polarizabilities 
in Table 1. 

(5) For a comparison of the ChPT - predictions of the dipole polarizabilities 
with calculations performed before 1994, and for additional information 
on these quantities, we refer the interested reader to Refs. [5] and [26]. 

7.3 Estimating the uncertainties 

To estimate the uncertainties in the prediction of the polarizabilities, we first 
note that the helicity non-flip amplitude H ++ starts out at order p 4 . We have 
therefore, for this quantity, a leading and next-to-leading order calculation at 
our disposal. For the corresponding polarizabilities (a\ — fti)-* and (a 2 — ^2)-^°, 
we thus simply add in quadrature the uncertainties generated by the order p 4 
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and p 6 LECs (see Section 3). The resulting numbers are given in column 2 of 
Table 1. They do not incorporate an estimate of the higher order contributions. 

On the other hand, the helicity flip amplitude H + _ starts out at order p 6 , and 
we have determined here only its leading order term. According to Eq. (2.7), 
this amplitude is proportional to B(s,t,u), which is an analytic function of 
the variables s, v at the Compton threshold and can be, therefore, expanded 
in a Taylor series, 



B(s, t,u) = U + Vs + Wu 2 + 0(s 2 , z/ 4 , su 2 ) . 
The relation to the polarizabilities is 



(7.7) 



(cti + Pi)^ = 8aM w U , (a 2 + fa)-** = QGaM^V . (7.8) 
The Taylor coefficients themselves have a chiral expansion of the form 



U = 



1 



V 



(16tt 2 F2) 2 
1 

(16tt 2 F 2 ) 2 M 2 



M 2 U t 
16tt 2 F 2 
M 2 V, 



V + 



16tt 2 F 2 



+ 0{M$) 
+ 0{Mt) 



(7.9) 



Whereas LECs from order p 6 do contribute to U , the leading term V is a 
pure loop effect, because Vq/M 2 is not analytic in the pion mass and thus 
cannot receive contributions from polynomial counterterms. To illustrate this 
point, and to estimate the size of V±, we consider the vector meson exchange 
amplitudes worked out in [5]. The contribution from uj exchange is dominant 
and given by 



C 

B^S.t.U) = -y 



+ 



Ml-t Ml-u 



; = 0.67 GeV" 



(7.10) 



In the language of ChPT , this amplitude starts out at order p 6 . It does con- 
tribute to Uo - this term is included in the resonance exchange estimates for 
the 0(p 6 ) LECs in (3.11). For this reason, we calculate the uncertainties in 
(«i + /3i) n o as before, with a result that is given in the second row of Table 1. 
Again, it doe not incorporate an estimate of higher order contributions. 

Finally, we come to the estimate of the uncertainty in (a^+Z^)?!- - m agreement 
with what is said above, resonance exchange does not contribute to Vq. On the 
other hand, there is no reason why this term should dominate the contribution 
from V\. Indeed, if we use (7.10) to also estimate effects from order p 8 , we find 
with M w = 782 MeV the value V\ = —2.2. The corresponding contribution 
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to the quadrupole polarizability (a 2 + /^vr is -0.25 • 10- 4 fm 5 - of the order 
needed to bring the ChPT calculation into agreement with the analysis of 
Ref. [12]. 

This result illustrates that the discrepancy between the chiral prediction for 
the quadrupole polarizability {0.2+ fa)-* at order p 6 and the dispersion analysis 
in Ref. [12] is of no significance, because the terms neglected may well be much 
larger than the leading order term, which would only dominate for very small 
values of the pion mass. On the other hand, to obtain a reliable estimate 
of (a 2 + P2)n° i n the framework of ChPT , one needs to perform a reliable 
calculation of the relevant couplings at order p s . This is outside the scope of 
the present work. For this reason, we do not quote an uncertainty for (a 2 +P2)-K° 
in Table 1. 



8 Summary and outlook 

(1) We have recalculated the two- loop expression for the amplitude 77 — > 
7T 7r° in the framework of chiral perturbation theory. We have made use 
of the techniques developed in Ref. [10], and of the effective Lagrangian 
Cq available now [8,9]. 

(2) The method has allowed us to evaluate the dipole and quadrupole polar- 
izabilities in closed form. [As far as we are aware, the quadrupole polar- 
izabilities have never been calculated in ChPT before.] The two Lorentz 
invariant amplitudes A and B are presented as a sum over multiple in- 
tegrals over Feynman parameters whose numerical evaluation poses no 
difficulty This is in contrast to Ref. [5], where part of the amplitudes, 
denoted by Aa,b, were published in numerical form only. 

(3) Our result agrees with the earlier calculation [5] up the the coefficient 
in one of the chiral logarithms in the amplitude A, and up to minute 
differences in the numerical values of the remainder Aa,b- The induced 
changes in the numerics of the cross section and of the dipole polariz- 
abilities are far below the uncertainties generated by the (not precisely 
known) values of the low energy constants. 

(4) The values for the dipole and quadrupole polarizabilities are presented in 
Table 1 and confronted with recent evaluations from data on 77 — > 7r°7r°. 
There is reasonable agreement for the dipole polarizabilities. As for the 
quadrupole ones, the combination (a 2 — fa)* related to the helicity non- 
flip amplitude agrees with [12] within the uncertainties quoted. On the 
other hand, the sum (a 2 + ^2)^° - related to the helicity flip amplitude - 
differs in sign from the one in Ref. [12]. We have shown why this does not 
contradict the predictions of ChPT : this quantity is a two- loop effect, 
and one expects from order p 8 (three loops) substantial corrections to 
the leading order result. We have indeed identified u;-exchange as an 
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important contribution at this order. 
(5) It would be instructive to improve the estimates for the LECs C29, . . . , C34 
in the sense that in these estimates, the constraints from the asymptotics 
of QCD [27,19] should be respected. 

The corresponding calculation of the charged pion polarizabilities is in progress 
[16]. 



Acknowledgements 



This work was completed while M.A.I, visited the University of Bern. It is 
a pleasure to thank G. Colangelo, J. Bijnens and B. Moussallam for useful 
discussions, and S. Bellucci for useful remarks concerning the manuscript. This 
work was supported by the Swiss National Science Foundation, by RTN, BBW- 
Contract No. 01.0357, and EC-Contract HPRN-CT2002-00311 (EURIDICE). 
M.A.I, also appreciates the partial support by the Russian Fund of Basic 
Research under Grant No. 04-02-17370. Also, partial support by the Academy 
of Finland, grant 54038, is acknowledged. 



A Notation 

In order to simplify the expressions, we set the pion mass equal to one in all 
Appendices, 

M n = l. (A.l) 

We use the following notation for <i-dimensional one-loop and two-loop inte- 
grals, 



<">=/<2?M->- «■■■»= /^i/^w- (A - 2) 

In particular, 



= F n [z],n>l, 



, [z ~ I 2 ] 



F n [z] = C(w) r(n r( ' } W) ,w = \ - 2 . (A.3) 
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The measures in the integration over Feynman parameters are defined by 



d 2 x = dx2dxs , d s x = dx\dx^dx^ . (A. 4) 

In dispersion relations, we use the <i-dimensional measure 



1 J T(3/2 + w) V4 / H ' 
^H=(4^' P= V 1 " 4 /- (A.5) 



B The acnode 



The evaluation of the two-loop vertex and box diagrams (5) and (6) in Fig. 3 is 
described in Ref. [10]. It is based on a rf-dimensional dispersive representation 
of the fish-type diagram, 



Jtf) = ( [1 _ P][1 1 _ {l + p) , ] ) = CM r<-«o jdx [i-fxd- x T 

= /J^L;-1.5<c<0, (B.l) 

where the measure [da] is given in Appendix A. This representation allows one 
to reduce the two-loop vertex and box integrals to the one-loop case, with a 
final integration over the dispersion parameter a. The ultraviolet divergences 
can be extracted by invoking recursion relations [10]. While the acnode was 
treated in a different manner in [10], we evaluate it here analogously to the 
vertex and box diagrams just mentioned. This results in considerable simpli- 
fications in the numerical programs. For this purpose, we invoke a dispersion 
relation for the function I(m, n; s) defined by 



I(m, n;s) = J dx [1 — s x (1 — x)] m [x (1 — x)] n , — 1 < m < . 
o 

(B.2) 

I(m,n;s) is analytic in the complex s-plane, cut along the real axis for Re 
s > 4. To evaluate its absorptive part, we observe that the imaginary part of 



18 



the first factor of the integrand in Eq. (B.2), evaluated at the upper rim of 
the cut, is 



Im [1 - sx (1 - x)] m = - sin(Trm) [sx(l-x)- l] m , 

s>4 x-<x<x+, x± = (1/2) ^1 ± . (B.3) 

The integrand is symmetric around x = 1/2, so we may restrict the integration 
from x_ to 1/2 in the evaluation of the absorptive part of I(m,n;s). The 

substitution x = (1/2) (l — ^(1 — 4/ s) (1 — tt)j generates a hypergeometric 
function, and we arrive at the dispersion relation 



1 (to, n; s) = / 

J a 

4 



da p (to, n; a) 



r(3/2)5 1 / 2+m /l 3 3 \ 

P(m ' n; a) = 4"r(3/2 + m)r(-m) 2 ^ U ' 2 + m + n; 2 + m] 'V ' 

5 = ^-1. (B.4) 
In particular, we find 



T(3/2)(3S r 

p(to,0;<7)- 



r(3/2 + TO)r(-TO) ' 

. 1 . T(3/2)f35 m / 2\ 

P(to, 1; a) = - — - — 1 + to + - , 

FV 7 4T(5/2 + TO)r(-TO) V a) 

. n . T(3/2)(35 m ( 2 4(1 + to) 12\ 

P(rn,2-a)= )'> M 2 + 3to + to 2 + ^ ) - + -\ , 

16 r(7/2 + m)T(— m) \ a a 1 ) 

(B.5) 

with p given by Eq. (A. 5). 

The dispersion relation (B.4) allows us to evaluate the integral that occurs in 
the evaluation of the acnode diagram (1) in Fig. 5, 



AflV 

N 



KK D 1 D 2 D z D 4 D 5 )h 



V L = (l 2 + q 2 -hY -1, 

V« = (/i + gi-/ 2 ) 2 -l, 

£>i = 1 - l\ , L> 2 = 1 - (Zi + qif , 
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D 3 = 1 - l\ , D A = 1 - (Z 2 + g 2 ) 2 , 
D 5 = 1 - (Zx - Z 2 - pi + gi) 2 . 



(B.6) 



We combine l/(-Di -D 2 ) and 1/(D 3 D4) by using :ri and x 2 as Feynman param- 
eters, respectively. By shifting Zi and Z 2 and again dropping terms that vanish 
upon contraction with the polarization vectors, one obtains 



1 



|/" - IcPrrU l l ll PN(k,Pi,qi) u 

r = h~Q, Q = Xiqi + x 2 q 2 - Pi ■ (B.7) 

Here, PN{hiPiiQ.i) is a polynomial in the momenta indicated. The integration 
over Zi is performed by using the dispersion relation (B.4) with s = r 2 . Then 
we proceed in a manner which is similar to the case of the box diagram de- 
scribed in Ref. [10]. The final expression can be written as a combination of 
the integrals 



A(i,k,m,n) = j[da] JcPx(^-l) * a~ k (1 - x 3 ) m F n [ Zacn ] , 
4 

z aC n = xl + (1 - x 3 ) cr - x 3 (1 -x 3 )a, 

a = Xi x 2 s + X\ (t — 1) + x 2 (u — 1) , (B.8) 

where F n [z] is defined in (A. 3). The integrals A(i,k,m, n) are convergent at 
w = in the case 



i = 1, fc = 0, m > 1, n > 4, 
i = l, fc > 1, m > 1, n>3, 
i = 2, A;>0, m > 0, n>3. 

To single out the divergent part in the remaining integrals, we invoke recursion 
relations in the following manner. We perform a partial integration in x 3 , and 
use 

[dx 3 (i-x 3 r = - {1 ~ x ^ +1 . 

J m + 1 

Then we express (1 — x 3 ) a through z acn and obtain 

(m + 3 — n + u;)A(i, fc, to, n) = 

Div(i, k,n)-n [A(i, k, m, n + 1) - (1 + a) fc, m + 2, n + 1)] , (B.9) 
where 
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Divf. k n)-e—^ g » r (3/2) Y{n-2-u) 
DlV(? ' /£ ' n) " 4 r(3/2 + c) r(n) 

5(5/2 -i + u;,n + &;-4-H- 2u) , 

The divergences in Div(i,k,n) can be worked out straightforwardly. 
The integral A(l, 0,0,3) must be considered separately. We write 

oo 1 

A(1,0,0,3) = L>(0,3) + |[Ax] J d 3 x{F 3 [ Zacn ] - F 3 [y]} , 

4 

y = xl + (l-X 3 )<T. (B.ll) 

The divergent quantity D(0,3) is worked out in Appendix C.l of Ref. [10], 
whereas the integral on the right-hand side is convergent at d — 4. 

This concludes our discussion of the acnode integral (B.7). 



C The quantities and A B 



Here we display the expressions for the quantities Aa(b) in Eqs. (5.2) and 
(5.3). 



A A (s,t,u) 
4 



689 4tt 2 \ 15043 



(4ttF w ) 4 [ ^162 9 J 64800 J 
{ F A n (s,t,u) + FT(8)+F%"( 8 ,t,u)}, 



(4ttF w ) 4 288 

1 / 8329 /2987 2vr 2 \ ll 
Ab ^ U > ~ j4^Fj* j 43200 + V 1350 " ~9~ ) ~s) 



(4ttF^) 4 288 



(C.l) 



(C.2) 



where 



pa.cn 



oo 1 



p(0) p(l) 
7; acn ^ 7; ac 



+ 



3(2) 
I; acn 



-acn ^ a cn 
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da 13 } 2 P /;vcr 



Z v 



'da (3 



Fi° X =J-^ J d 3 - E {P^L + D { ^ + + PSox_££L } ; / = A, B , 

4 n=1 

(C.3) 



and 



£>W - 1 ± 1 

bOX± 7T7u * 7 n 

^"<box;t ^box;u 

Here Pj are polynomials in s, z/ = t — u and in Xj. Their explicit expressions 
are given in Appendix D. The quantity z acn is displayed in Eq. (B.8), y is given 
in (B.ll), and 



z ver = a (1 - x 3 ) + xly 2 , y 2 = l-sx 2 (l- x 2 ), 

z box; t = B t — A t X\, 

A t = x 2 x 3 [s(l-x 2 )x 3 + (1 -t) (1 -x 3 )], 
Bt = A i + z ver , 

^box; u ^box;t|t— *u • (^"^v 

The acnode integrals are easy to evaluate numerically in the physical region 
for the reaction 77 — > 7T7r, because branch points occur at t = 4, -u = 4 only. 

On the other hand, the vertex and box integrals contain branch points at 
s — 4. In order to evaluate these integrals at s > 4, we invoke dispersion 
relations in the manner described in [10]. 



D The polynomials P4 and P B 



Here, we display the polynomials Pa{b) that occur in the expressions Aa(b) in 
Appendix C. We use the abbreviations 



x + = x\ + x 2 — 2 x\ x 2 , x_ = x\ — x 2 , 
2723 = (1 + a; 3 - 2x2X3) (1 - x 3 + 2xix 2 x 3 ) . (D.l) 
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D. 1 The polynomials Pa 



PfLn = ~ 192^3 (1 - Xg)( S X + - VX-), 



p aLu = -6svx-(1-x 3 ) 



X 



1 + 8 x 3 + 3 x\ - 4 x + (x 3 + 3 x 3 3 + 2 xj) + 2 x 2 _ (1 + 2 xl - 3 xf) 
6 s 2 (1 - x 3 ) [x^ (1 + 2x 3 + 14x 3 + 7x1) +6x 2 _xl 
1 + 2 x 2 . + 4 ^ (2 + x 2 _) + 3 x\ (1 + 2x! 



-x 



+6 z/ 2 xl (1 - x 3 ) 1 - 2 x 3 + 2 X3 + 5 Xg - 12 x+Xg 
+12 s [-2 x^ (1 - x 3 ) 3 (1 + 2x 3 + 3x^) 

+Xg (2 + 33 Xg - 19 Xg - 19 X^ + 15 Xg) 

+ X + (-4 + 2 X 3 - 63 Xg + 25 X^ + 51 Xg - 35 Xg) 

+2 xt x 3 (2 + 9 x 3 - 7 Xg + 5 x\ - 3 x 3 ) 
+12z/x_ (1 -x 3 ) [4-2x3-5x3-4x3 + 19x3 



21 x^; 



19 x\ + 10 x A 3 



+2x + (I + 8X3 
-48 x 3 4 + 4x 3 + 6x3 
+x+ (2 - 3 x 3 - 19 X3 + 41 Xg - 21 x 3 ) 

^A; ) acn = 2sZ/X_(l-X 3 ) 2 [-11 -8x 3 



8^3 



24 x 3 



9 x 3 



+ 12 x+ (4 + x 3 + 3x 3 i + 2 (11 + 8x3 + 8x3 + 6x 3 - 9x3) 

-2 s 2 (1 - x 3 ) 2 [x^(35 - 18 x 3 + 34 x 3 + 21 x\) + 18 x 2 _ x 4 3 
/ ll + 16x 3 + 9x^ + 2x 2 _(ll 



-x. 



2x 3 



$4: 



+2 v 1 xt (1 - x 3 y [-13 + 18 x 3 - 2 x^ + (15 - 36 x+) 
+4 s (1 - x 3 ) [x 3 (22 + 11 x 3 + 83 x 2 - 101 x\ + 45 x£) 
-2x 2 + (l-x 3 ) 2 (ll(l+x 3 + x 2 ) -9x 3 ) 



.1 o 



+x + (60 - 178 x 3 + 43xg 



173 x^ + 233 xl 



105 x^) 



+2x 2 _(2-x 3 )x 3 (ll(l + x 3 + x 2 ) + 9x: 



+4z/x_ (1 - x 3 ) ; 



56 x^ + 57x^ 



-60 + 74x 3 + 9x 3 : 
+x+ (22 + 104 x 3 - 126 x^)" 

+ 16 x 3 (1 - x 3 ) [60 + 18 x 3 - 108 x 2 3 + 93 x^ - 30 x A 3 
+x+ (-22 - 83 x 3 + 205 x 2 3 - 187 x 3 + 63 Xg) 
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P A; vcr = 128 sx\ (1 - 2 x 2 ) x\ 3 - 2 x 2 (2 - 15x 3 )x 3 
-32s 2 x 2 (l - 2 x 2 ) (30 - 42 x 2 + 7X 2 ,) x^ , 



18 x 2 



15 x 2 



A; boxa 



-96 s x 2 x 2 [—6 + 9 x 3 + 20 x 2 x 3 — 22 xi x 2 x 3 



—2 x 2 — 40 x 2 x 2 + 6 x\ x 2 x\ — 4 x 2 x 2 + 52 xi x\ x\ — 8 x 2 x 2 x 2 



-15 Xn 



X 2 Xg 



14 Xi x 2 x 



12 



VO rp ,7-* 2 , T ,3 o 2 2 3 



+8 xi x 2 X3 + 4 x : x 2 X3 + 20 x 3 — 48 x 2 x 3 — 36 x\ x 2 x 3 



+76 xi x^ X3 + 10 x? x 2 X3 - 24 xi x^ Xg - 24 x{ x\ x\ + 24 x{ x\ x\ 
-16 s 2 x\ Xg [l2 + 12 xi - 12 x 2 - 24 x x x 2 + 12 x\ x 2 + 30 x x x 3 
—24 x 2 x 3 — 90 xi x 2 x 3 + 27 x 2 x 2 x 3 + 30 x 2 x 3 + 66 X\ x\ x 3 

QO o o o o 

—24 x : x 2 x 3 — 10 x'l x 2 x 3 — 24 x 3 + 6 Xi x 3 + 66 x 2 x 3 — 18 X\ x 2 x 3 



„2 ^,3 



20 x\ x\ 

„2 _4 _4 



2 2 
c l •^2 1 

„2 „ J2 



+45 x^ x 2 Xg - 48 xi, Xg + 48 x l x\ Xg - 120 x\ x\ Xg - 28 x\ x 2 x 2 
-24 xi + 36 

*^ 1 ^2 ^3 

+ 56 

*^ 1 "^2 "^3 



+48 z/ 2 x^ Xg 



—4 + 8 Xi — 8 x 2 + 2 xi x 3 + 7 x 2 x 3 + 10 x 2 x 3 



—30 xi x 2 x 3 + 30 Xi x 2 x 3 — 10 x'l x 2 x 3 + 6 x 3 — 14 x x Xg + 3 x 2 x 2 

— 12 X 2 Xg + 12 Xl X 2 Xg + 12 X 2 X 2 Xg — 12 X^ X 2 Xg + 24 Xl X 2 Xg 

AS -I- 94 1"^ 

*Aj Y *L 2 ^g I ZjT: ^1 **> 2 *^ 3 

+384 x 2 x 2 —6 + 9 x 3 + 19 x 2 x 3 — 19 xi x 2 x 3 + 7 x 2 — 50 x 2 x 2 

O 0*3 *3 '-I 

— 2 Xl X 2 Xg + 52 Xl X 2 Xg — 15 Xg + 44 X 2 Xg + 16 Xl X 2 Xg 



60 Xl X 2 Xg + 8 Xg — 16 X 2 Xg 



16 xi x 2 Xg + 32 xi 

3^2 ^3 



p(2) 

-'A; box + 



—24 s v 2 (1 — Xi) 2 x\ (1 + x 2 + Xi x 2 — 2 Xi x 2 ) Xg Xi 23 
96 s x 2 Xg x i23 — 7 x 2 — 7 xi x 2 + 14 xi x 2 + 3 x 3 + 6 x 2 x 3 



+6 xi x 2 x 3 — 12 xi x 2 x 3 — 3 x 3 — 2 x 2 x 3 



2 xi x 2 x 2 + 4xi ^£3 



—48 s x 2 x 3 xi 23 — 2 — 2 xi — x 2 + 6 xi x 2 — x x x 2 + 3 x 3 + 3 xi x 3 
—6 xi x 2 x 3 — 4 x 2 — 4 xi x 2 + 2 x 2 x 2 + 8 xi x 2 x 2 

, n 1 2 a 2 2 /i222_|_/i23 2 

I £1 X-j^X 2 Xg 4r Xl X 2 X g TzX-j^X^Xg | 4r X ^ X 2 Xg 

-24 s 3 x 2 (1 + xi - 2 xi x 2 ) (1 + xi - x 2 — x 2 x 2 ) Xg x i23 

-48 v 2 (1 - xi) 2 x\ (1 - x 3 ) x\ (1 + 4 x 3 ) xi 23 + 2 x 2 (1 - x 3 ) 2 Xg 1 xi 23 



A; box_ 



325Z^3J2^3 



— 6 + 6 xi + 12 xi x 2 — 18 Xi x 2 
— 15 xi x 3 + 12 x 2 x 3 + 33 x 2 x 2 x 3 — 30 xi x 2 x 3 + 21 

^ 1 ^2 ^ 3 

— 10 x 3 x 2 x 3 + 12 Xg + 3 xi x 2 — 24 x 2 x 2 — 30 Xi x 2 x 2 

+ 18 x 2 x 2 x 2 + 6 x 2 x 2 + 54 xi x 2 x 2 — 48 x 2 x 2 x 2 + 4 x\ x\ x\ 
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-96 v x\ x\ [-9 + 11 xi + 12 x 3 - 14 x 1 x 3 + 24 x 2 x 3 
—56 Xi x 2 x 3 + 28 x\ x 2 x 3 + 15 x 2 — 27 Xi x 2 — 58 x 2 x 2 
+84 x\ x 2 x\ + 8 x 2 ^2 X3 + 60 x\ x\ x\ — 80 x\ x\ x\ 
-16x1 + 28 xi xl + 32x 2 xl - 20x1X2 X3 - 42x 2 x 2 x\ 
-64x1X2X3 + 84x^X2X3 , 

P A;box_ = 48 S V I 1 - ^l) X \ A ^123 

x — 2 + 3 x 3 + 3 X2 X3 + 3 xi X2 X3 — 6 xi X2 x 3 — 4 X3 

—2 x 2 x 2 — 2 xi x 2 x 2 + 4 xi x\ X3 

+24 s 2 v (1 - Xi) X2 (2 - x 2 - xi x 2 ) X3 X123 
+u (1 - xi) x\ (1 - x 3 ) X3 (1 - 2 x 3 ) X123 

+24 z/ 3 (l -Xi) 3 X^ X%X X23 . 



D.2 The polynomials Pb 



B; acn 



96 x 3 (1 - x 3 ) 



x_ 



X+ 



p(l) 

B; acn 

3z/ 2 



61s 



X- 



4 - 2 x 3 + 21 xi - 19 x^ + 3 X3 + 5 x 



s 
24 



xl(l-x 3 ) [l + 2x 3 + 2x 3 + x|] 



-— x 3 (2-x 2 ) [2-8x 3 + 5x^] +3s(l - x 3 ) 
s 

x x\{l- x 3 ) 3 (1 + x 3 ) - x + (1 + 8 x\ + 3 4) +6i 2 _i^ 
+6 [x + (1 - x 3 ) 2 (-4 - 2 x 3 - 3 x 2 + 13 xl) 
+x 3 (-2 - 33 x 3 + 19 x 2 + 19 x 3 - 15 xj) 



>(2) 



2v 



B;acn--— X- (1 - X 3 ) 

-34 + 152 x 3 - 43 xl + 21 x^ - 27 x 4 3 + 15 x^ 



x 2 . (1 - x 3 ) 2 [35 - 18 x 3 - 2 x 3 + 3 x 4 3 

+-x 3 (l -x 3 ) [34 - 67x 3 + 95x 2 - 56x 3 + 15x 3 t 
s 

-s (1 - x 3 ) 2 [x 2 + (1 - x 3 ) 2 (13 + 8x 3 + 3x^) 
+x+ (11 + 16 x 3 + 9 x*) - 18 x 2 x 4 3 
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+v x_ (1 - x 3 ) 2 11 + 16 x 3 3 + 9 xl + 12 x + (4 - 3 x 3 - x|) 
+2(l-x 3 ) [-x 3 (22 + 11x3 + 83^ - 101^ + 45 x\) 
+x + (34 - 108 x 3 + 65x 2 3 + 73 x\ - 103 x 4 + 39 x\) 



P B; vcr = 16sx 2 (l - 2 x 2 ) x^ [33 x 2 x 3 - 6 (4 - 5 x 3 ) - 2 x 2 (4 + 9x 3 ) 
-1536 x^l - 2 x 2 ) (1 - x 3 ) (1-2 x 3 ) x\ , 



P, 



(i) 



24 z/ 2 



B; box+ 

1152 



x^ (1 — x 3 ) X3 4 + 8 x\ — 10 x 3 — 2 xi x 3 — 3 x\ x 3 



i 2 (1-13)13 (I-213 + 24) 



-24 s x 2 xj [4 + 12 xi - 4 x 2 - 12 x 3 - 14 x 1 x 3 + 10 x 2 x 3 

—14 x\ x 2 x 3 + 17 x\ x 2 x 3 — 14 x\ x\ x 3 + 8 x\ — 2 x\ x\ — 6 x 2 x\ 



+ 14 xi x 2 X3 — 3 x\ x 2 X3 + 4 xi X2 x 2 — 12 x^ X2 X3 + 12 x\ x\ x 3 



„2 2 2 



„2 3 2 



48 x 2 (1 — x 3 ) x 2 —6 + 19 X3 + 16 x\ x 2 X3 



27x1 



2 x 2 x\ 



— 18 Xi x 2 x 2 — 28 x\ x\ x\ + 20 X3 — 32 x\ x 2 x\ + 64 x\ x\ x\ 



p(2) 

- r B;box + — 
1 . 

+ -x 2 (l 



24 z/ 2 



(l-x 1 ) 2 x3(l-x 3 ) 3 x|(l + 4x 3 ) 



3 ^4 



x 3 ) 4 x^ 



/ \ O /I o 

— 24sx 2 (l — X3) x 3 — 2 — 2 xi — x 2 + 6 xi x 2 — x x x 2 

+3 X3 + 3 xi X3 — 6 xi x 2 X3 — 4 x 2 — 4 xi x 2 + 2 x 2 x 3 

+8 Xi x 2 x 2 + 2 x 2 x 2 x 2 — 4 x\ x\ x 3 — 4 x 2 x 2 , x 2 + 4 x\ x\ x\ 



-12 X2 (1 + xi - 2 xi x 2 ) (1 + xi - x 2 - xj x 2 ) (1 - x 3 ) 2 x 6 
+ 12 z/ 2 (1 - xi) 2 x\ (-1 - x 2 - xi x 2 + 2 xi x 2 ) (1 - x 3 ) 2 X3 
+48x 2 (l -x 3 fx\ [-7x 2 



7 Xi x 2 + 14 xi x 2 + 3 x 3 

2 ™ q _2 o ™ ~2 



+6 x 2 x 3 + 6 Xi x 2 x 3 — 12 x\ x 2 x 3 — 3 x 3 — 2 x 2 x 3 
—2 xi x 2 x 2 + 4 xi x 2 X3 



P (i) 

-'B; box_ 



48 v 



x 2 , (1 - x 3 ) 2 X3 (5 + 11 xi - 28 x 3 - 4 xi x 3 ) 



+48 v x 2 x 3 [2 + 2 xi + 8 x x x 2 — 6 x 3 — 3 X\ x 3 — 2 x 2 x 3 



—4 xi x 2 x 3 + 3 x x x 2 x 3 — 12 X\ x 2 x 3 — 3 x 1 x 2 x 3 + 4 x 3 

It 



+Xi x 2 + 2 x 2 X3 — 10 x 1 x 2 x 2 



2 2 

1 37]^ 



p(2) 
-'B; box_ 



144 v 



(l-x 1 )x 2 2 (l-x 3 Yxt(l -2x 3 ) 
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+—(l-x 1 ) 3 x 4 2 (l-x 3 ) 2 xl 

-12 s v (1 - x\) x\ (-2 + x 2 + xi x 2 ) (1 - x 3 ) 2 

+24 z/ (1 — xi) x 2 (1 — X3) 2 X3 [—2 + 3 X3 + 3 X2 X3 + 3 xi X2 X3 

—6 Xi x\ x 3 — 4 x 2 — 2 x 2 x 2 — 2 xi x 2 x 2 + 4 xi x 2 X3 . 
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